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Abstract 

We study the transverse momentum resummation for dijet correlation in hadron collisions based 
on the Collins-Soper-Sterman formalism. The complete one-loop calculations are carried out in 
the collinear factorization framework for the differential cross sections at low imbalance transverse 
momentum between the two jets. Important cross checks are performed to demonstrate that the soft 
divergences cancelled out between different diagrams, and in particular, those associated with final 
state jets. The leading and sub-leading logarithms are identified. All order resummation is derived 
following the transverse momentum dependent factorization at this order. Its phenomenological 
applications are also presented. 
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I. INTRODUCTION 


Dijet production in hadronic collisions is one of the golden channels to study perturbative 
QCD and hadron physics in high energy experiments. In dijet events, the two jets are 
produced mainly in the back-to-back conhguration in the transverse plane, 


A -\- B —y Jct\ -1“ J 6^2 “1“ X 


( 1 ) 


where A and B represent the two incoming hadrons with momenta P and P, respectively, the 
azimuthal angle between the two jets is defined as 0 = 0i — 02 with 0 i ^2 being the azimuthal 
angles of the two jets. There have been comprehensive analyses of the azimuthal angular 
correlation (or decorrelation) in dijet events produced at hadron colliders In the 

leading order naive parton picture, the Born diagram contributes to a Delta function at 0 = 
TT. One gluon radiation will lead to a singular distribution around 0 = vr, which will persist at 
even higher orders [^. This divergence arises when the total transverse momentum of dijet 
(imbalance) is much smaller than the individual jet momentum, = l-Piu + Aul ^ l-Pi±l ~ 
l^2±| Pt, where large logarithms appear in every order of perturbative calculations. These 
large logs are normally referred as the Sudakov logarithms, a* ln^*“^(P^/g^). Therefore, a 
QCD resummation has to be included in order to have a reliable theoretical prediction. 
The goal of this paper is to derive the resummation formulas for dijet production in the 
kinematics of these large logarithms. A brief summary of our results has been published in 
Ref. @. 

In the kinematics region q± Pt, the appropriate resummation method to apply is 
the so-called transverse momentum dependent (TMD) resummation or the Collins-Soper- 
Sterman (CSS) resummation [^. The CSS resummation was derived for the massive neutral 
particle production in hadronic collisions, such as the electroweak boson (or Higgs boson) 
production. However, due to the presence of the colored final state, the resummation of 
dijet production will be much more complicated. There have been theoretical progressed 
in the leading double logarithmic approximation (LLA), where it was found that each in¬ 
coming parton contributes half of its color charge to Sudakov resummation factor [M 


In this paper, we will go beyond the LLA to perform the resummation calculation at the 
next-to-leading logarithm (NLL) level. The threshold resummation for dijet production in 
hadronic collisions has been investigated in a series of publications by Sterman and his col¬ 
laborators [ll|, 12|. The methodology of our calculations for the TMD resummation is very 


similar to those studies. 

We will start our derivation by evaluating the differential cross sections as a function of 
the imbalance transverse momentum q± at the complete one-loop order. The leading order 
contribution is a Delta function 5^‘^\q±). One-loop corrections come from four contributions: 
(a) virtual contributions; (b) soft gluon radiation (real); (c) jet contributions; (d) collinear 
gluon radiation associated with the incoming parton distributions. The virtual graphs have 
been studied in the literature [l^. The jet contributions are easy to derive following the 
examples of inclusive jet production. In order to calculate the analytical results, we will 
adopt the narrow jet approximation (NJA) jH, 14|, where explicit dependence on jet sizes 
Ri and R 2 can be evaluated. The soft gluon (real) contribution is most difficult to calculate. 
This is because, not all soft gluon radiation contributes to the finite q±. The radiation inside 
the jet will be part of the jet contribution at one-loop order and has to be excluded from soft 
gluon radiation contribution. In our calculations, we will apply a small offshellness for the 
jets in the hnal state to calculate the soft gluon radiation. We hnd that this method yield 
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the exact same results as that using kinematic cut-off to regulate soft gluon radiation in the 
NJA. Collinear gluon radiation associated with incoming parton distributions also contribute 
to the hnite imbalance transverse momentum. This part can be formulated according to the 
well-known DGLAP splitting. In the end, for hnite and soft imbalance transverse momentum 
q±, we add the soft gluon (real) and collinear gluon (real) contributions together, which leads 
to the so-called asymptotic behavior for the differential cross sections at low imbalance 
transverse momentum q±. An important cross check of our derivation is to compare its 
numerical results with the dijet production codes available in public. We will carry out 
these comparisons in this paper. Another cross check is to demonstrate the soft divergences 
cancellation between real and virtual graphs. We will show these cancellations in details for 
the hard partonic channels in our calculations. After the cancellation, we are left with only 
collinear divergences associated incoming parton distributions. This indicates that we do 
have a consistent results for dijet production at one-loop calculations. 

The large logarithms mentioned above will be evident from the complete one-loop results. 
Resummation of these large logarithms is the main goal of current paper. To perform 
the resummation, we hrst show that the differential cross sections at one-loop order can 
be factorized into the transverse momentum distributions, soft factor, and hard factors, 
respectively. The transverse momentum distributions follow the dehnitions for Drell-Yan or 
Higgs boson production at low transverse momentum. The soft factor will have to take into 
account the additional effect of gluon radiation associated with the two hnal state jets. The 
idea to construct the soft factor follows the examples of threshold resummation studied by 


Sterman et ah 11 


Our resummation formula can be summarized as 




dyidy2dP^(Pq_i 


5^ CTO 

ab 


d% 


(27r 


M r 


Wab^cd{Xl,X2,hA_) + Yc 
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where the hrst term W contains all order resummation and the second term Y comes from 
the hxed order corrections. (Tq represents the overall normalization of the differential cross 
section, yi and y 2 are rapidities of the two jets, Pt the leading jet transverse momentum, 
and gx fhe imbalance transverse momentum between the two jets as dehned above. All 
order resummation for W from each partonic channel ah —)■ cd can be written as 


Wah^cd {Xi,X2,h) 


= Xi fa{xi, /i = ha/h2_)x2 fb{x2, jJ. = &o/&±)e ^Sud(Qhfex) 


X Tr 


'_A 

Safe^cdexp 

rQ 



Jbo/b± h 


. dfeo/fex 


da , 

—7 

/i 


( 3 ) 


where = s = X 1 X 2 S, representing the hard momentum scale. 60 = 2 e , with 7 ^ being 
the Euler constant. fafi{x,y) are parton distributions for the incoming partons a and b, 
and Xi ^2 = Pt (e^^^ -f e^^^) /y/S are momentum fractions of the incoming hadrons carried 
by the partons. In the above equation, the hard and soft factors H and S are expressed as 
matrices in the color space of partonic channel ab —?• cd, and 7 ^ are the associated anomalous 
dimensions for the soft factor (dehned below). The Sudakov form factor Ssud resums the 
leading double logarithms and the universal sub-leading logarithms. 


^sud(Q^&±) = 




d/i 
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where Ri ^2 represent the cone sizes for the two jets, respectively. Here the parameters A, B, 
Di, D 2 can be expanded perturbatively in as- At one-loop order, A = B = —2Ca/5o^ 

for gluon-gluon initial state, A = B = ^ for quark-quark initial state, and A = 

(Cf+Ca) ^^ q ^ _ CaM^- for gluon-quark initial state. Here, /^o = (H - 2Ar^/3)/12, 

with Nf being the number of effective light quarks. At the next-to-leading logarithmic 
level, the jet cone size enters as well j^. That is the reason we have two additional factors 
in Eq. (|1]): D 12 = Ca^ for gluon jet and D 12 = for quark jet. The cone size R is 

introduced to regulate the collinear gluon radiation associated with the final state jets. Only 
the soft gluon radiation outside the jet cone will contribute to the imbalance q± between the 
two jets. 



It was found, in particular, certain diagrams in dijet production in hadronic collisions can 
not be factorized into the simple universal TMDs. In terms of resummation coefficients, this 
will affect the coefficient in Eq. (jl]). In our numeric calculations, we will include both 
A^^'> and coefficients in the resummation formula (see detailed discussions in Sec. HC). 
Though we do not expect that the factorization violating effect will affect much the results 
to be presented below, it will be useful to estimate the size of factorization breaking effect 
in some future work. 

Second, we derive our results, in both collinear and TMD approaches, by adopting the 
narrow jet approximation. With that, we are able to demonstrate the explicit cancellation of 
soft divergences from various parts of analytical calculation, and to express the resummation 
Sudakov form factor in an analytic form in terms of the jet sizes. Although the derivation 
itself is not limited to NJA, we hnd that the differential cross section expressions are much 
simplihed under the NJA. How to extend our results to the case without NJA is very 
interesting subject. It is worthwhile to come back to this question in the future. However, 
in current paper, we will carry out our calculations with the NJA. 

The rest of this paper is organized as follows. We will dedicate Secs. H-V to the calcu¬ 
lations of the differential cross sections at low q± in the collinear factorization approach. In 
Sec. H, we start with a brief introduction of the leading order results in which the overall 
normalization factor for each production channel is presented. In Sec. HI, we discuss the 
generic features to evaluate the one-loop corrections to Wipp). This includes virtual graph 
contributions, jet contributions in the real gluon radiation, collinear gluon radiation asso¬ 
ciated with the incoming partons, and the soft gluon radiation which contributes to hnite 
g_L. Since the soft gluon radiation is the most important contributions for dijet production, 
we dedicate Sec. IV to show the details of our calculations. In Sec. V, we compare the 
asymptotic results for dijet production to the fixed order calculation published in the litera¬ 
ture. The asymptotic behavior for dijet production is calculated by summing the soft gluon 
radiation and collinear gluon radiation derived in Secs. HI and IV. In Sec. VI, we derive the 
TMD factorization and resummation. In particular, we will show that the collinear calcula¬ 
tions from Secs. HI-IV at the one-loop order can be factorized into the TMDs, hard and soft 
factors. The latter factors are expressed in a matrix form, where they are written in terms 
of the color space for every partonic channel. The resummation is achieved by solving the 
relevant renormalization equations. In Sec. VH, we perform the phenomenological studies 
based on our resummation formalism. We conclude our paper in Sec. VHI. 
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FIG. 1. Schematic diagram for dijet production at the leading order, with incoming parton mo¬ 
menta: Pi and P 2 -, and outgoing momenta: ki and A: 2 , respectively. 

II. DUET PRODUCTION AT THE LEADING ORDER 


Dijet production at the leading order can be calculated from partonic 2 —)■ 2 processes, 

a{pi) + b{p2) -)■ c{ki) + d{k2) , (5) 

where pi ^2 and ki ^2 are momenta for incoming and outgoing two partons. Schematically, we 
draw the diagrams in Fig. [T] 

The partonic channels include the following subprocesses: 


qtqj qiqj , 
qiqi qjqj , 
qiqi -t qiqi , 

qg ^qg , 
gg^qq, 
qq^ gg, 
gg ^ gg ■ 

Their contributions can be summarized as 
d^a 


dyidy2dP^d‘^q± 


2 ^ 2 . = ^ 03^1 fa{xi,p)x2 fb{x2, {qj 

ab 


,( 0 ) 


( 6 ) 

( 7 ) 

( 8 ) 
( 9 ) 

( 10 ) 

( 11 ) 

( 12 ) 

(13) 


+ 


where the overall normalization of the differential cross section is The partonic 

cross sections for all the production channels are listed below. 

(14) 

(15) 

( 16 ) 

( 17 ) 

(18) 

(19) 

( 20 ) 
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FIG. 2. Schematic diagrams for virtual graph contribution (a) and final state jet contributions (b) 
and (c) at one-loop order. Both of them are proportional to a Delta function of the imbalance 
transverse momentum: S^‘^\q±). 


where the kinematic variables s = (pi+P 2 )‘^, t = {pi — ki)‘^ and u = (pi —^ 2 )^. As mentioned 
above, at the leading order, they contribute to a Delta function of q±, which corresponds to 
the back-to-back conhguration of the two jets in the transverse plane. If we translate this 
into the fo^-space, we will obtain W{h_[_) at the leading order take the following form, 

= Xifaixi, p)x2fb{x2, p)h^aLcd ' (21) 

The goal of the following three sections is to derive the one-loop corrections to 


III. GENERIC DISCUSSIONS ON ONE-LOOP CALCULATIONS 

The leading order contributions lead to a Delta function of the imbalance transverse 
momentum In the following, we will first carry out one-loop calculations of the differential 
cross sections at low transverse momentum. The computation is performed in the fe^-space, 
the Fourier transformed conjugate parameter to the transverse momentum q±. That is to 
say, we will calculate W{h±) at the one-loop order. 

There are four types of radiative contributions at the one-loop order: 

1. Virtual contributions, as shown in Fig. [2](a). These contributions have been calculated 
in the literature, and they are proportional to the leading order contributions, leading 
to a Delta function of g_L. 

2. Real gluon radiation: Jet contributions. In the real gluon radiation, one particular 
contribution is also proportional to the leading order contribution, that is the jet 
contribution. In this case, the gluon radiation is within the jet, where its momentum 
is collinear to the hnal state parton. For example, as shown Fig. [2](b), the radiated 
gluon is collinear to one of the hnal state parton, and they form a new jet at one-loop 
order. Because of momentum conservation, this again leads to a Delta function of q±. 

3. Real gluon radiation: collinear gluon associated with incoming partons. These 
collinear gluons contribute to a hnite g_L, and is proportional to l/q\ multiplied 
by the splitting kernel of parton distribution functions. Since they originate from the 
incident partons, the collinear gluon contributions follow the generic structure and are 
easy to calculate. 

4. Real gluon radiation: soft gluon contribution. Soft gluon contributions are more 
difficult to evaluate. They also contribute to a hnite q±. To evaluate this part of 
contribution, we apply the leading power expansion in the limit of g_L Q. However, 
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FIG. 3. Schematic diagrams for real gluon radiation contribution to finite imbalance transverse 
momentum: (a) collinear gluon radiation associated with the incoming partons; (b) soft gluon 
radiation outside the jet cone of final state jets. 


because the final state jets also carry color, the soft gluon radiation has to take into 
account the interactions among initial and hnal state partons. In addition, we have to 
exclude the soft gluon radiation within the hnal state (cone) jets whose contributions 
have already been included in the hnal state jet contributions. Detailed calculations 
will be presented in the following section. 

In the rest of this section, we will go through the hrst three kinds of contributions, whereas 
the soft gluon radiation contribution will be calculated in Sec. IV. 

We will carry out our calculations in the collinear factorization formalism, and apply the 
dimensional regulation for IR and UV divergences with dimension D = 4 — 2e. Various 
divergences appear in individual contributions: 1/e^ represents the soft divergence, whereas 
1/e for either soft or collinear divergence. Since we are dealing with jet production in the 
hnal states, we will also encounter the divergences associated with the jet size R. The 
explicit calculations of gluon radiation will depend on how we dehne the jet, i.e., the jet 
algorithm will play a role in formulating the one-loop corrections. 

Two important cross checks will be performed in the derivations. First, the soft diver¬ 
gences of 1/e^ will be cancelled out completely among diherent terms. We notice that (1,2,4) 
terms listed above will have 1/e^ contributions. A crucial test of our calculations is that 
these 1/e^ cancel out each other. This cancellation is not trivial, since they come from dif¬ 
ferent diagrams, and some with diherent color factors. However, we will show that the total 
contribution is free of soft divergence of 1/e^. Second, the divergences associated with the 
hnal state jets are also cancelled out among various terms. The collinear divergences asso¬ 
ciated with the jets are regulated by the jet sizes, and the individual contributions contain 
terms of l/eln(l/i?) are cancelled out in the hnal results. We will show that the cancellation 
indeed happens at this order. 


A. Virtual Graphs 

Virtual diagrams for dijet production in hadronic collisions have been calculated in the 
classic paper of Ref. [l^ . In our calculations, we will take their results. 


B. Jet Contributions 

Jet contributions contain the collinear gluon radiation and gluon to quark-antiquark 
splitting in the hnal state. The schematic diagrams have been shown in Figs. [2](b) and 
(c). Because we have two jets in the hnal state, one-loop jet corrections can come from 
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either of the jets, as shown in Fig. |2]^b) and (c), respectively. The jet contribution needs to 
be included to capture the collinear gluon radiation (or gluon splitting to quark-antiquark 
pair) within the jet cone. These gluon radiations will not change the kinematics of the 
parenting parton, and therefore contribute to a Delta function of q±, which is similar to 
the virtual graph contributions. The requirement is that the two partons in the splitting 
process form a jet according to the jet algorithm adopted in experimental measurements. 
In order to derive an analytic expression with the jet cone size dependence to demonstrate 
the cancellations in the hnal results and to derive the resummation formula, we apply the 
narrow jet approximation (NJA) in our calculations. In particular, we follow the technique 


and scheme of the subtraction developed in Refs. [13|, . The basic idea is to note that in 


the collinear gluon radiation limit for each hnal state jet. 


\M{2 3)p ^ \Mo{2 2)|2 x V 


1^-2 


( 22 ) 


where Ad(2 —)■ 3) represents the scattering amplitude of 2 —)■ 3 subprocess, M.q{2 2) for 

the leading order 2 —)• 2 subprocess with one of the hnal state parton branching into two 
parton hnal state represented by the splitting of Vi^ 2 - 
Therefore, the jet contributions can be summarized as. 


f d^h Ej^ 


(23) 


where the factor Ej/E 2 accounts for the phase space diherence from 3 parton hnal state to 
2 parton hnal state. The phase space integral of the above equation is limited that the two 
partons are within the jet cone. Here the diherence from jet algorithms plays a role. Using 


the NJA, the calculations follow what have been done in Refs, 
we hud that the gluon jet in the hnal state contributes 


and in particular. 




Os 1 1 

^r(l -e) -2e 


PtR\ 

dvr / 





UggiO + fqg{0] , 


(24) 


for fcrtype jet algorithm, where R dehnes the jet cone size: R = ■ Here, A?/ 

and A0 are the rapidity diherence and azimuthal angle diherence between the two partons 
which dehne the jet. In the above equation, fgai£) and fqg{^) are splitting kernels for gluon 


to gluon and gluon to quark-antiquark pair [R 
subtraction, we obtain. 


1^ . respectively. By applying the MS 


J^ = 


ttgC/ 

2ti 


1 1 , P^R^\ 


J— 2/?o ~ In ■ 




) 


+ p 


(25) 


where we have taken into account the contributions from both g ^ gg and g ^ qq splittings. 
We would like to emphasize that the singular terms (double and single poles) are independent 
of jet algorithm. The jet algorithm dependent contributions arise from the hnite term P. 
The above jet function is universal, and can be used in many other cases as well. Similarly, 
we hnd that the quark jet can be written as 




agCp 

271 




/i 2 ) 


+ P 


(26) 
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FIG. 4. Example of Feynman diagrams that contribute to the TMD factorization breaking effects 
in two particle production via hadronic processes at 0{al) order [l^ (left), where two gluon 
exchanges between the incoming/out going particles induce non-universality of the TMDs in this 
process as compared to other processes such as DIS and Drell-Yan processes. These diagrams 
would contribute to in the resummation formula. The dominant contribution at this order 
is illustrated in the right panel, which are factorizable and come from in the resummation 
formula. 


In the /crtype jet algorithms, the above mentioned and P terms are already available: 


P 

P 


In —h— 
P 

In —r— 
/i 2 


- 2/?o In 

P 

-In- 

2 /r 2 


P^P2 gy 3 23 ,, 

+ — - -vr^ - —N; 
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54 


/ ’ 


+ 


13 


4^ 


(27) 

(28) 


We note that the double logarithmic terms are independent of jet algorithm, and it is 
important for them to cancel in the hnal results, as shown in Sec. IV. 


C. Collinear Gluon Radiation 


The contributions from the collinear gluon associated with the incoming partons can be 
easily evaluated, and they are found to be proportional to the splitting kernel at the one-loop 
order: 


0 (c 


dx[ dx2 


27r2g2 


x[ 


X 2 


fa'{x[)fb'{x2) + (^1 ^ 2 ) 


/a 


(29) 


where = Xijx^ and P^^j/ represents the splitting kernel part without the Delta function 
term whose effect is included in the above mentioned virtual contributions. 

To evaluate the contribution from soft gluon radiation requires more care. This is because 
we have to exclude the collinear contributions associated with the hnal state jet which will 
be factorized into a jet function and will not contribute to hnite q±. Before we discuss its 
details in the following section, we note that the TMD factorization breaking effects can 
appear at higher orders in two particle production at hadron colliders [i3-l2il. These effects 
come from the diagrams illustrated in Fig. S] [l^, which belong to a nontrivial contribution 
at order to the parton distributions. They cannot be factorized into the conventional 
transverse momentum dependent parton distributions, though they could contribute to the 
coefficient in the resummation formula. Since it is beyond the perturbative order (up 
to a^) discussed in this paper, we shall not discuss it further in this work. 
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FIG. 5. Soft gluon radiation contribution to finite imbalance transverse momentum q±: (a) initial 
state radiation, and (b) and (c) final state radiations. Because we have chosen the gluon polarization 
vector along p 2 , there is no contribution connecting to p 2 line. 


IV. SOFT GLUON RADIATION AT ONE-LOOP ORDER 


For soft gluon radiations, we can apply the leading power expansion and derive the dom¬ 
inant contribution by the Eikonal approximation. This analysis has been applied in Ref. 
to obtain the leading double logarithmic contributions to dijet production. In the current 
paper, we will extend that analysis to include the subleading logarithmic contributions as 
well. In particular, we will derive the contributions which depend on the jet cone sizes. The 
relevant Feynman rules have been listed in Ref. [8| . For completeness, we copy these results 
here. For outgoing quark, antiquark and gluon lines, we have 

2 fcf 

2ki ■ kg + ie^ ’ 2ki ■ kg + ie^ ’ 2ki ■ kg + ’ 

respectively, where fcj represents the momentum of the outgoing particles. For incoming 
quark, antiquark and gluon lines, we have, 

2pi „ 2,p1 2p‘^ 


-9 


-9 


-9 


(31) 


2 pi ■ kg — ie 2pi ■ kg — ie 2pi ■ kg — ze‘ 

respectively, wher^i represents the momentum for the incoming particle. 

Following Ref. |8|, we choose the physical polarization for the soft gluon along the incom¬ 
ing particle p 2 - Therefore, the soft gluon radiation from the incoming particle p 2 vanishes 
with this polarization choice. From that, we can derive the soft gluon radiation contribution 
easily, with the polarization tensor for the radiated gluon. 


^^^{kg) = -9^^ + 


k^P^2 + 

kg ■P2 


(32) 


fll/ 


flL' 


For example, from the amplitude squared of the soft gluon radiation terms in the above, we 
have 

2p)^ 2p'i . 

2pi ■ kg 2pi ■ kg 

2k^ 2k'i . 

2ki ■ kg 2ki ■ kg 

2k^ ‘^k^2 ^ 

2k2 ■ kg 2k2 ■ kg 

where S'g(p, q) is a short-handed notation for 

Sg{p,q) = - 


5'g(Pl,P2) , 

(33) 

Sg{kl,P2) , 

(34) 

Sg{k2,P2) , 

(35) 

■ 9 

1 ■ ko ' 

(36) 
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Similarly, we derive the interferences between them, 

^Pl 


2kl ■ kg 2pi ■ kg 
2k2 ■ kg 2pi ■ kg 

2k^ 2k% , 

2ki • kg 2k2 • kg 


fii/ 


flip 


Sg{kl,P2) Sg{pi,P2) - Sg{kl,Pl) , 

(37) 

Sg{k2,P2) + Sg{pi,P2) - Sg{k2,Pl) , 

(38) 

Sg{kl,P2) + Sg{k2,P2) - Sg{ki,k2) . 

(39) 


In order to evaluate the contributions from soft gluon radiation, we integrate out the phase 
space of the gluon whose transverse momentum leads to the imbalance between the two jets, 

e.g., 

o f d^-^kn 


9 


-6^^\q± - kg±)Sg{pi,P2) , 


(40) 


. {27r)D-^2Ek^ 

where we have chosen dimensional regulation for the phase space integral. The derivation 
of the above term is straightforward, by noticing that the lower limit in the longitudinal 
momentum fraction integral, 

dx 1 


/ 'r 

^ ^min gJ- 


(41) 


where we have dehned x as momentum fraction of pi carried by the soft gluon. Because of 
momentum conservation, we have lower limit for the x-integral: Xmin = Therefore, the 
above integral leads to the following leading contribution. 


1 , 

91 91 


(42) 


Substituting the above equation into Eq. fl40|) . we have 
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2 


d^ ^kg 

{2n)D-^2E,^ 


d^^\9. 


kgx)Sg{pi,P2) 


Ois 1 
2 'k‘^ q\ 



(43) 


Because there is no e = (4 — Zl)/2 term in the dx integral, we have taken H = 4 in the 
above equation. The other terms are more difficult to calculate, because the above phase 
space integral will contain jet contributions which have already been taken into account by 
the above mentioned jet functions. To avoid double-counting, we have to subtract this part 
of contribution. That means the phase space integral has to exclude the jet (cone) region. 
It is interesting to note that the amount of this exclusion does not depend on jet algorithm. 
This is because, here we are considering the soft gluon radiation, whereas the jet algorithm 
mainly focuses on the collinear gluons associated with the jet. 


A. Out of the Jet-cone Radiation 

As a general discussion, let us take the example of one term. 
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where ki± represents the transverse momentum for the hnal state jet, and = kg-pi/ki -pi. 
Clearly, there is collinear divergence associated with the jet. That means, if the gluon 
radiation is within the jet cone, it will generate a collinear divergence. In order to regulate 
this collinear jet divergence, we can limit the phase space integral to require that the gluon 
radiation being outside of the jet cone. With this restriction, there will be no divergence 
associated with the jet. Instead, the jet (cone) size R will be introduced to regulate the 
collinear divergence from the jet. 

There are different ways to regulate the above integral. The main task is to identify the 
need of introducing jet (cone) size. In the above example, the integral diverges when kg is 
parallel to ki, where the invariant mass of ki + kg becomes small. The out of cone radiation 
requires that the invariant mass has a minimum, say. A, i.e., 

(fci + kgf > A' . (45) 

Clearly, A depends on the jet size. In other words, if {ki + kg)"^ is smaller than A^, we have 
to exclude its contribution, because it belongs to the jet contribution calculated in previous 
section. 

Following the similar analysis as done for the jet contribution, we can hnd out the size 
of A. For example, if we substitute the kinematics of fci and kg into the above equation, we 
will obtain 


{ki + kgf ^ ki±kg± - 2ki±kg± cos(0i - (j)g) Ri ki±kg±Rlg , (46) 

where yi and Pg are rapidities for ki and kg, 0i and (pg are the azimuthal angles, respectively, 
and Rig represents the cone size between ki and kg. In other words, if Rig is smaller than 
R, the gluon radiation will be considered inside the jet cone. Therefore, in the phase space 
integral of Eq. fl4T)) . we have to impose the following kinematic restriction: 0(2fci ■ kg — A^) 
with h? = ki±_kg±R^. Equivalently, we hnd it is much easier to adapt a slight off-shell-ness 
for the jet momentum ki to regulate the divergence: k\ = nil — kl^R?. By doing that we do 
not need to impose any kinematic constraints, and the phase space integral is much easier 
to carry out. The choice of ml is to make sure that {ki + kgY is always larger than hi. This 
can be verihed as follows, 

{ki + kgY = \lklj_ + mlkg_i_ - 2ki±kg± cos(0i - (pg) 

ijklj_ + ml - kiA^ . (47) 

By choosing ml = kl^R^, it is guaranteed that {ki + kg^ is larger than A^ for any values of 
/hy and hep. 

In the narrow jet approximation, i.e., in the i? —?• 0 limit, the phase space cut-off technique 
results into the same leading contributions in terms of ln(l/ R). After adding an off-shell-ness 
to the jet momentum, the integral in Eq. fl4T)) can be written as 


\Jkl^ + mlkg^{/hyY + fci±fcg±(A0)^ -F 2kg± 


j ^kg±S^^\q± 



^_ 2 _ 

6 eiklA^ + R^) + - 2iiki^ ■ V ■ 


(48) 


To proceed, we average the azimuthal angle of the jet but £x the azimuthal angle of kg±. 
This corresponds to keeping the imbalance transverse momentum direction q±. With that. 
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we obtain 


r d4 _ 2 _ 

J Jo TT + R^) + k‘^g_^-2^iki±kg±cos{(j)) ' 

fd^ _ 1 _ 


( 49 ) 


where the lower limit of is an illustration, we have taken D = 4 in the average 

of 0 angle in the above equation. There is a e-term correction if we keep D-dimension, which 
can be calculated accordingly. In the hnal results shown below, we have kept those terms 
for completeness. After taking the limit of Q and i? —)■ 0, we obtain the leading power 
contribution as 



In ^ + In 
Q1 




(50) 


Therefore, the hnal result for the integration of the Sg(ki,pi) term can be written as 


g 


^^3 kgA_)Sg{ki,pi) = 


{‘2nf2Eu 


27r2 ql [ ql 


In ^ -h In ^ -h In ( - 
u 


Rl 


Evaluation of the other terms are similar, we summarize their hnal results as 

ds 1 


Sg{.PuP2) 

Sg{ki,Pi) 

Sg{k2,Pl) 

Sg{kl,P2) 

Sg{k2,P2) 

Sg{ki,k2) 


271 ^ ql 


21 n^ 

Ql 


27r2g2 

1 

l-t^l 

ds 1 


27r2g2 

Hi Q 

L Ql 

1 

In"^' 

27r2g2 

Hi Q 

L Ql 

ds 1 

L, 

27r2g2 

1 

H 

1-to 

ds 1 

r, 1 


In —^ 


Rf 


In —h In —77 + In I — ) + € ( — In —77 H—— 


Rl 


u 


u 


TT 


Ri 6 


In —; 5 —h In —77 -|- In ( e I — In -\—— 


TT 


Rl 


271 ^ ql 


Rl 


u 


In —-|- In —77 -|- In ( — ) -|- e - In —77 -|- — 


TT 


2 Rl Q 


Rl 

1 


In ^ -h In --77 -h In - + e { - In -h — 


u 


TT 


2 ■■■ Rj 6 


Rl 


+ In ^ + 2 In 




Rl 


\tu 


+ 6 


1 , 2 ! 1,2 

- In^ ^ -I— In^ 

2 Rl 2 


Rl 


+ 


7T 


4 1 -5 1 -S 

-4 In — In- 

—t —u 


(51) 

(52) 

(53) 

(54) 

(55) 

( 56 ) 

(57) 


where we have kept the e terms which will make hnite contributions to the complete one- 
loop calculation of W{h), when Fourier transforming from to fe^-space. Furthermore, the 
leading double logarithm contributions arise from the A- In ^ terms in the above equations. 

The above results are the basic elements to be used in our calculation for deriving the 
low g_L behavior of a scattering process, induced by soft gluon radiation. In the following, 
we will apply these results to all the partonic processes which contribute to the inclusive 
dijet production in hadronic collisions. 
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B. qiqj qiqj 

Quark-quark channel (with different quark flavors i and j) is the simplest case to calculate. 
Its leading Born amplitude can be written as 

Mo = u(/ci)T“7^u(Pi)u(/c2)T“7''u(P2) {G^y{ki - pQ) , (58) 

where represents the gluon propagator with momentum ki — pi. For simplicity, we 
separate the color factor in the above amplitude, 

Mo = Ao{pi,p2, ki, /c2)hfcT“uiU/T“Mj , (59) 


where i,j and k, I are color indices for the incoming and outgoing quarks, respectively. Soft 
gluon radiation can be summarized as 


2n^ 

- ^AquT^T^uuT^u 

2pi ■ kg 

2k^ 

+ ^, \ AoumuT^u 

2ki * kg 

2k^ 

+—-^AouT‘^uurT^u , 
2k2 • kg 


(60) 

(61) 

(62) 


where kg is the radiated gluon momentum, p, and c for its polarization vector and color 
index, respectively. To calculate the soft gluon contribution via this partonic channel, we 
need to perform the phase space integration over its amplitude squared, as discussed in the 
previous subsection. 

Let us hrst work out the color factors for various terms in its amplitude squared: 


k>tk'i^GF\Mo\^ , 
k^k]^ ^ Gf\Mo\^ , 
pW ^ Gf\M,\^ , 

Kkl => -\Gf\M,\^ , 

ki^p\^-^{2GA-GF)\Mo\^ 


(63) 


Including the proper color factors, the amplitude squared contributes 


Gf [Sg{Pl,P2) + Sg{ki,P2) + Sg{k2,P2)] 

+ 0^ + Sg{pi,P2) - Sg{kl,Pl)\ 

Zl\c 

~ {‘2Ga - Gf) [Sg{k2,P2) + Sg{puP2) - ^ 3 (^ 2 ,pi)] 

-^GF[Sg{k2,P2) + Sg{k,,P2) - Sg{k,,k2)] - (64) 
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After integrating the (restricted) phase space of the radiated gluon, we obtain the following 
contribution of the soft gluon radiation to the qiQj —)■ qiqj channel; 


cts 1 
27r2 q\ 




n2 

2Cf + 
T± 




-^2/ J 



(65) 


where 


F (<?<?') 
sn 




2 {Ca - Ci.) In — - 2 In — 
—t —u 


( 66 ) 


An important cross check of the above result is to show that the infrared divergences of the 
soft gluon radiation are cancelled by those from the virtual diagrams and jet contributions. 
The only left divergences are associated with the collinear divergences for the incoming 
two quark distributions, which can be absorbed into the dehnition of renormalized parton 
distribution functions. 

To check the cancellation, we have to Fourier transform the above expression into the 
impact parameter fe^-space. 


^ I {ftg 


X 


X 


Xi X 2 
1 &0 

1 bl 

■ 7 ^ ""Tv 


qiqj^qtqj 


Cf ( -«') + K « o) + ^(i - 0^(1 - a 


Cpln 


p2 p2 

K 1 K 2 


+Cf 




fI 2 ^ 1 , 2 ! 

+ ‘'^^( 2 '" b! + 2‘“ + T 

—t —u 


In" 


\ bl J 6 

+5(1 - 0^(1 


( 67 ) 


where we have also included the collinear gluon radiation contributions associated with the 
incoming two quarks. 

The virtual graphs have been calculated in the literature, and can be summarized as 
follows. 


CXg 


Cf 


4 1 

FT H- 






u 


4 In —— 2 In —^ + 2 In-6 

/i^ t 


1 1 ,, s" 

H-4 In — 

e 2Nc tu 


+ ••• 


( 68 ) 


where we only kept the singular terms to check the cancellations between real and virtual 
diagrams. In addition, we have two jets contributions 


CK, 


Jeti + Jet 2 = —Cf 
27r 


— H— I 3 
e 


p2 

2 In ^ + In 


p2 p2 

iXliX2 


+ 


(69) 


where P* are hnite terms associated with jet functions. Clearly, the 1/e^ and 1/e terms all 
cancel out after summing up all the above three contributions, except those associated with 
the splitting of quark distribution: 



<^(1 “ i') + ('C + 35(1 - 05(1 - 0) 


= -1 iF,(?)i(i - ?')+- «)i. 


(70) 
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where Vqq is the quark splitting kernel. The complete expression for the hnite terms will 
be much involved. To facilitate the discussion on factorization, to be presented in Sec. VI, 
we show below the most important terms in the hnite contributions, in particular, those 
with large logarithms of and ln(/i^65_/^o)- This will clear show how the TMD 

factorization works. 




- hi (p„(mi - (')+ 


(f) 


+ In 


6S 


—?>Cf + Ci? In —2 + Cp In 

Ki K 2 


-5(1 - 05(1 - 0) In 


Q^bi 

bi 




(71) 


where an overall integrand factor of J fq{x[)fq{x2) was omitted for simplicity. We 

would like to emphasize a number of important observations from the above calculations. 
First, the factorization scale fi dependence only exists in terms associated with the parton 
splitting kernel. This scale dependence shall be cancelled by the relevant scale evolution 
for the integrated parton distributions. Second, in the hnal results, ln^(l/i?^) terms are 
cancelled out between the jet contribution and the soft gluon contribution. Third, the large 
logarithms appear in the one-loop calculations contain three terms: (a) the double logarithms 
in terms oi\n‘^{Q‘^b‘]_/bl) proportional to incoming partons color factors (here, it is Cp + Cp); 
(b) single logarithms in terms of Inlfp^b^/bl) associated with parton distributions; (c) the left 
single logarithms of ln{Q‘^b\/bl) contains similar terms as Drell-Yan process (the —3Cp term) 
and those associated with dijet production in this particular channel (jet size dependent 
contributions and additional contributions which is process-dependent). All these features 
point to a possible factorization in terms of TMDs, for which we will discuss in Sec. VI. 


C. qg qg 

In this process, we have two different color structure at the Born level, 

AiuT^T\ + A2 uT^T^u , (72) 

where a and b represent the color indexes for the incoming and outgoing gluons, the ampli¬ 
tudes Ai and A 2 depend on momenta of two incoming particles: pi and p 2 , and two outgoing 
particles: ki and k 2 for the quark and gluons, respectively. The leading order amplitude 
squared reads as, 

|Ao|' = Cp (Ai + A 2 )' - CaAiA; , (73) 

where the two terms are separately gauge invariant. Soft gluon radiation follows previous 
example, and can be decomposed into the following three terms, 

[AimT"T“T'’u + A2urT^T^u] 

i-tfcbd) [AiuT'^T'^u + A2uT‘^n] , (74) 


2k^ 

2ki ■ kg 

-2K 

2pi ■ kg 

2kli 

H- - — 

2k2 ■ kg 
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from the initial state and final state radiations, where c represents the color index for the 
radiated gluon. 

The amplitude squared of the soft gluon radiation can be summarized into the following 
form, 


[CFSg{pi,P2) + CFSg{h,P2) + CASg{k2,P2)] 


^ (^1 + 


-^{Ai + A2f + -2AiAl 


If. . .2 a2 

- {A,+ A2f + ^Al 


{Sg{kl,P2) + Sg{k2,P2) - Sg{ki,k2)) 

{Sg{kl,P2) + Sg{pi,P2) - Sg{kl,Pl)) , 
{Sg{k2,P2) + Sg{pi,P2) - Sg{k2,Pl)) • 


(75) 


Adding them together and applying the phase space integral, we obtain the leading contri¬ 
bution induced by soft gluon radiation in the qg —)■ qg channel; 


J_ J /i(o) 

27r2 ql I 


(Ca + Cf) ln^ + CF\n^ + CA\n^ 

-^1 -^2 J 


+ , 


where represents additional contribution in the sub-leading logarithm. 


rto) = In- 

—u 

+ In ■ 


2{N^ + l)s^ + u{s^ + u‘^) 


m 


—su 


+ 


—t‘^s 


-Nr, 


-I- 
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+ 


-su 




(Ca - Cf) 


(76) 


(77) 


The Fourier transformation of the above results into the impact parameter fe^-space leads 
to the following contributions. 


V,,•(.+.) (jn = ^ J ,,(o) 


27r 




-1 + In 2T1 (CFV„(mi - o + CAV„(e)s(.i - o) 
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In" 
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-7 + 

Virtual graphs contribute to the following terms in hF(&_L), 


g .1 ftW) 


1 1 
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In — 
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Furthermore, the jet contribution, including both quark and gluon jets in the hnal state, 
yields 

W'Hb.) = g [(C. + C.) (i - 1 ,n + i (|c. + + C. in ^ + C. in ^)|o) 

Clearly, all the divergences are cancelled out between the above terms, except the collinear 
divergences associated with the incoming quark and gluon distribution functions. 

Again, the finite contributions take the following form, if we only keep the logarithmic 
terms. 




Cf + Ca ^^2 
2 


- e ')+- 0] - <5(1 - o<5(i - r) 


- 2/3o + In ^2 + C'a In ^2 


Q% 1 \ ri°»' 

*0 / /iSt,, 


D- 99 99 

Similarly, the Born amplitude for the gg —)■ qq channel is 

AiuT^T% + A2uT^T% , (82) 

with two momenta for incoming gluons: pi and p 2 , and two momenta for outgoing quark 
and antiquark: ki and k 2 . The leading order amplitude squared can be written as 

= Cp + A 2 ) — CaAiA2 , (83) 

with crossing symmetry to the above qg —)■ qg channel. Soft gluon radiation can be derived 
as 

2k^ 

-^ + A2hT‘=T^T“nl 

2ki - kg'- 

—2k^ 

+ -^ 

2k2 * kg 

2r)^ 

+ i-^fcad) [A.uT^^T^ + A2UT^T‘^v] , (84) 

Zpi • fVg 

where c represents the color index for the radiated gluon. Its amplitude squared, including 
proper color factors, yields 

|AoP [CASg{pi,P2) + CFSg{kl,P2) + CFSg{k2,P2)] 

~ (Al + ^ 2 )^ + - 2 A 1 A 2 {Sg{kl,P2) + Sg{k2,P2) - Sg{kl,k2)) 

■ 1 1 

+-+ -(Al + ^ 2 )^ {Sg{kl,P2) + Sg{pi,P2) — Sg{kl,Pl)) , 

■/V2 1 1 

- ^A 2 --(Ai + A2)^ {Sg{k2,P2) + Sg{pi,P2) - Sg{k2,Pl)) , (85) 
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After applying the integral over the phase space of the radiated gluon, we obtain 
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( 86 ) 


(87) 


The above result shows the leading contribution at low imbalance transverse momentum q±. 
After Fourier transformation, we have the following contribution to hF(6_L): 


hF(*+'=)(6^) = 
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( 88 ) 


The virtual graph contribution can be summarized as 




( 0 ) 

gg^qq 


-2Cp-2CA){-^--\n%)+- 
' e /iV e 


Nr In 


—tfu 2u/ 
t 


+ In ■ 


-u f t 


u 


2R 

TT 


+ 


-3Cf — 2/3o) 
1 + 


tu 


In 


tu 


.(89) 


Again, we hnd out the soft divergences cancel between the above terms and the jet contri¬ 
bution which is the same as that of Eq. ([69]) • Here, the remaining collinear divergences are 
associated with the incoming gluon distributions. 

Keeping only the logarithmic terms, we obtain its finite contribution to W^^\b±) as 
follows. 
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E. gg gg 


For gg —)■ gg channel, we can simply write down the following decomposition for the Born 
amplitnde, 

■^ifabefcde -^2facefbde F ^sfadefb ce •) (91) 

where a, b, c, d are color indices for the glnons, with momenta for incoming glnons: pi and 
P 2 , and for ontgoing glnons: ki and ^ 2 - The leading order amplitnde sqnared can be written 
as 


|Ao|" = {Al + + AiA; - AiA; + A2A;) 

One soft glnon radiation takes the form, 


(92) 
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2Pi 

T 7 ) 1~ fgo.f \.Al ffbefcde T A2ffcefbde T A^fjdefbce\ 

^Pl ■ ^g 

The amplitnde sqnared of the above radiation can be written as, 
\Ao\^Ca [Sg{pi,p2) + Sg{kl,P2) + Sg{k2,P2)\ 

+ {S,{k„p^) + S,{h,p^)-S,{h,h)) ^ 


(93) 
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After integrating the phase space, we obtain the leading contribntion for soft glnon radiation 
in this channel as 
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where 
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(96) 


Again, to check the cancellation between different contribntions, we perform the Fourier 
transformation to the impact parameter-ft^ space, and hnd 
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where we have also included the collinear gluon contributions, and Vgg{^) is t he g luon-gluon 
splitting kernel. The virtual graphs have been calculated in the literature 1^, and they 
contribute, 
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The two outgoing gluon jets contribute 
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Comparing the above results, we hnd that the divergences do cancel each other except the 
collinear divergences associated with the gluon splitting functions. 

If we only keep the logarithmic terms, we will hnd the hnite contributions as 
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V. ASYMPTOTIC BEHAVIOR AND COMPARE TO FIXED ORDER COMPU¬ 
TATIONS 


An important cross check we will carry out in this section is to compare the asymptotic 
result of the dijet differential cross section in low region to the hxed order computa¬ 
tion derived in the literature. Through that, we could be sure that we have captured the 
most important contributions in the low q± limit. As discussed in Sec. Ill, the dominant 
contributions in the low q± region arise from collinear gluon radiation associated with the 
incoming partons and soft gluon radiation discussed in Sec. IV. Hence, the asymptotic result 
is obtained by adding these two contributions together. 

From the derivations of Secs. Ill and IV, we hnd that at low transverse momentum the 
diherential cross section has the following generic form: 
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( 101 ) 


where Ca, Cb, Cc, and Cd are the associated color factors for the incoming and outgoing 
partons: U = Up for quark and C = Ca for gluon. The calculations in the last section 
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FIG. 6. The comparisons between the asymptotic derivations of Eq. (jl01|) and the full LO calcu¬ 
lations for the kinematics specified by the DO collaboration at the Tevatron. The LO curves come 
from Ref. [l|. 


have presented the result of for some partonic channels. All other channels can be found 
accordingly. 

In Fig. ini we compare the above differential cross section, for dijet production at the 
Tevatron for the kinematics specihed by the DO Collaboration [ij, to the hxed order pertur¬ 
bative calculation with one gluon radiation contributions. From these plots, we can see that 
the above asymptotic results capture the leading contributions at low transverse momentum 
where 0 is close to tt for the back-to-back azimuthal correlation region. We would like to 
emphasize that the jet size dependent terms are crucial to make these comparisons. Without 
them, the agreements between the full LO calculations and our asymptotic derivation results 
will not be as evident as shown in Fig. |6l 

VI. FACTORIZATION AND RESUMMATION 

Clearly, from the derivations in the last section, and the plots of the differential cross 
sections shown in Fig. El we see that the collinear factorization calculations for the dijet 
production lead to divergent behavior at low imbalance transverse momentum q±, which 
corresponds to the back-to-back azimuthal angular correlation between the two hnal state 
jets. In this region, to have a reliable prediction, we need to perform the resummation. 
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Since we are interested in the dependence of the differential cross sections on the transverse 
momentum q±, the TMD resummation is the appropriate framework for carrying out our 
theory calculations. 

The TMD resummation was originally derived for Drell-Yan type of hard processes, where 
the hnal state only contains color neutral particle with large invariant mass. In order to apply 
this resummation formalism to the current case of dijet production in hadronic collisions, 
some modihcations are needed. In particular, the two hnal state jets carry color so that 
additional soft gluon radiation will introduce large logarithms associated with the hnal state 
jets. Although, these additional interactions will not modify the leading double logarithmic 
contribution, they will enter at the next-to-leading logarithmic level. This is evident from 
the explicit calculations discussed in the previous section. 

To carry out the TMD factorization, we extend the original CSS formalism, and take 
into account the hnal state radiation by assuming a soft factor in the factorization formal¬ 
ism. The exact same idea has been applied to resum large logarithms associated with the 
threshold logarithms for dijet productions. The similarity between the TMD and thresh¬ 
old resummation is not surprising, because they both deal with soft and collinear gluon 
radiations. 

In this section, we will argue the TMD factorization for dijet production based on the 
explicit one-loop calculations presented in the previous sections. The factorization is verihed 
in the transverse momentum space and 6_L-space at the one-loop order. At this order, we have 
to take into account the matrix form of the hard and soft factors. The explicit comparisons 
support the factorization we argued in this paper. 

With the explicit form of the TMD factorization formulas, we derive the resummation 
result by solving the relevant renormalization group equations. In particular, the TMDs 
follow the examples of color neutral particle productions, which have studied extensively in 
the literature. The additional soft factor obeys the renormalization group equation controlled 
by the anomalous dimension. We calculate the soft factor at one-loop order, from which 
we derive the anomalous dimension at one-loop order. Solving this renormalization group 
equation, we resume the sub-leading logarithms in dijet production processes. 

Generically, the low transverse momentum q± originates from collinear gluon radiation 
associated with the incoming partons and the soft gluon radiation effects. In particular, 
the collinear gluon radiation from the incoming partons contribute to the total transverse 
momentum of the hnal state particles. This is the dominant contribution in dijet production 
in hadronic collisions. The collinear gluon radiation associated with the two hnal state jets 
will not contribute to the total transverse momentum of the dijet. This is because they are 
factorized into the jet contribution which do not contribute a hnite q± as we have shown in 
Sec. III. On the other hand, soft gluon radiation among the incoming partons and hnal state 
jets will contribute to the low imbalance transverse momentum of the dijet. To account for 
this part of contribution, we introduce the soft factor in the TMD factorization. Since the 
hnal states carry color, the soft gluon radiation is more complicated than that for the color 
neutral particle production in hadronic collisions. 

We follow the method developed for threshold resummation in dijet productions in 
hadronic collisions in Ref. [ll|, where the soft factor is formulated in the orthogonal color 
space and both the incoming and outgoing partons are represented by Eikonal gauge links. 
Because of light-cone singularity in the parton distributions, we choose oh-light-front gauge 
links for the incoming partons. Meanwhile, the Wilson lines associated with the hnal state 
jets are constructed in such a way that only out-of-cone gluon radiation contribute to a 
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nonzero q±. Consequently, the soft factor in our calculations will depend on the jet cone 
size 0 . 

Following the above argument, we can write down the factorization formula as, 
df^CT f 

dyidy2dP^cPq_ -E ^0 J d k\^d k2_\_d Xafk2±_) 


ab 


xTr [Hab^c(i(Q^)Safe^cd(A±)] S‘'‘^\ki± + k2± + A_l — q±) , ( 102 ) 

where fa{,Xi,ki±) and fb{,X2,k2±) are TMDs and will be introduced in the following. Here, 
the hard factor and soft factor are expressed in the matrix forms in the color spaces for the 
incoming and outgoing partons. We can also express W{b±) in the fc^-space as 

Wab^cd = Xi faixi, b±, Cl, /i^ P^)X2 fbix2, Ci, 

X Tr [B.ab^cdiQ'^, P^ yi - 2/2, -Ri, R2)Sab^cd{b±, 2/1 - 2/2, -Ri, -R2)] ( 103 ) 


where we have shown all the explicit dependence of the TMDs and hard and soft factors. 

In the following, we will hrst introduce the TMDs, and then formulate the soft factors 
for all partonic channels. With these factors calculated in perturbation theory, we will show 
the above factorization is valid at one-loop order, by comparing to the derivations in the 
last few sections. 


A. Transverse Momentum Dependent Parton Distributions 


In the factorization formula, xfi{x, b±, C^, fR) is the Fourier transformation of quark TMD 
parton distribution xfi{x,k±X‘^-,X‘)- We follow the Ji-Ma-Yuan scheme [2^ to dehne the 
TMDs. For the quark distribution, we have 0 


Xfq{x, Px,C^/i^) = 


di-d'^U - 

-e 


eP+^ +ifcx-CL 


X 


F+(27r)3 

(P|vl/(C-, a) 4 (r, C±) 7 +^.( 0 , 0 a)d/( 0 )|P> 


{01 ( 6_l ; 00) (cx); 6 _l ) (0; cx)) (00; 0) 10) 


(104) 


For gluon one, it can be written as: 




X 


p|ri(r.a)cL(rXx)7+c.i»=(o,Oi)f*-+(o)|p, 

r r (105) 

{0|£l^ft,(6_L; oo)£))fe,^(cx); 6±)£„afe(0; (X))C.^bc{^\0 )|0) 

In the above equations, the relevant gauge links have to apply, in the fundamental and 
adjoint representations for the quark and gluon distributions, respectively. The gauge link 
is chosen along the direction n, 

-0 


, C±) = -Pexp —ig / d\v ■ A{\v -h C) 


(106) 


^ It is possible to factorize the soft factor in our paper into a soft factor and a jet function for the final state 
jet. By doing that, the soft factor may not depend on the jet cone size. We leave this for a future study. 
^ Here, we follow the original definition in Ji-Ma-Yuan 22|, where the soft factor is subtracted from the naive 


gauge invariant TMDs. If other subtraction method would be used, the associated soft factor definition 
would have been changed as well. 
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for both cases, in the appropriate representations of SU (3), adjoint for the gluon distribution 
and fundamental for the quark distribution, respectively. Similarly, we dehne the TMDs from 
hadron P, which depend on the gauge link along the direction v. The vectors v and v are 
off-light-front: v = (n"'', 0j_) with v~ 3> v~^ and v = (h’*', 0j_) with v~ -C h’*', and 

. Since non-light-like vectors exist, the TMD parton distribution depends on a new large 
scale = (2n • P)/v‘^ and a free parameter = ^^-1+ ■ Based on these dehnitions, we can 
calculate them under perturbative QCD order by order, and can be expressed in terms of 
the integrated parton distributions. At one-loop order, the gluon distribution can be written 


/g(x,6,C,/i,p) = 




s'-'A 


271 


dx' ,, 
—- < In 


X' 


^ [PggiOfgix'^p) + Vgq{Ofgix',p)] + - l)fg{x,p) 


X 




(inp^ -F 2/3o - l) 




1 , 2^ C!_£!_ 7 

2 " bi "^2 "^2 2 2 


(107) 


where ^ = xjx'. For the quark distribution. 


fq{x,b,C,f^,p) = 


agCp 1“ dx' 
271 


X' 


^ \PggiOfg{x\ /i) + Vqq{i)fq{x', /i)] + - l)/g(x, p) 


X 


In ) In 


bl 


- - In' 




71 

“T “ ^ 


In the above equations, Vij are splitting kernels. At one-loop order, we have 


Vgg{x) = 6 


1 — X 


+ 


X 


T^gq{x) = 


X (1 — a;)- 
4 [ 1 + {1 — xY' 


-1-a:(l - x)/9o5(l - x) ) , 


X 


^ / X 4 / 1 -|- x^ 3 ^ . i 

= 3 + 2^(1 - A ) . 


(108) 

(109) 

( 110 ) 
( 111 ) 
( 112 ) 


Similarly, we can obtain the fc^-dependent expressions for the TMDs from these references. 


B. Color Space Decompositions for the Soft and Hard Factors 


It has been shown in Refs. |ll|, for soft gluon radiation in dijet production, it is more 


convenient to construct the factorization in the color space matrix, where the soft and 
hard factors can be calculated in the orthogonal color bases. In particular, the associated 
anomalous dimension for the soft factors can be formulated and the relevant resummation 
can be performed accordingly. 

In this paper, we consider dijet production through the partonic 2 —)■ 2 subprocesses. 
Therefore, for each partonic channel, we will construct the color-space bases depending 
on the color indexes of two incoming partons and two outgoing partons. The soft gluon 
radiation does not modify the fundamental scattering structure for each partonic channel, 
so that the color-space bases are constructed for all order perturbative calculations. In 
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addition, the color-space bases are not unique. In our calculations, we follow those used in 


Ref. [ill]. For the quark-quark scattering subprocesses, we have fundamental representation 


from both incoming and outgoing partons, for example. 


i i i + td, 


€ + 91 ^ + i‘d, € + « ^ 9^ + ii 


^3 






(113) 


where j and k indicate the flavors of the quarks, and a, b, c and d are the color indices. For 
this channel, we have two independent color conhgurations. 


^ab r j: ^ab rpa'rpa' 

^Icd — ^acObd 5 ^2cd — -^ac-^bd ’ 


(114) 


corresponding to the color-singlet and the color-octet couplings, respectively. Similarly, we 
have the same decomposition for identical quark-quark scattering subprocess, 

qa + <ib^<ic+qd, qi + Qb^qi + ql- (ns) 


For gg —)■ qq channel. 


9a + 9b ^ qc + qd^ 

we have three independent color bases, 

^ab rpa' rpa' ^ab jabc'rpc' ^ab • £abc' rjic' 

^Icd — -^ac-^cd 5 ^2cd — ^ cd •> ^Zcd — cd ' 

Similarly for qq gg channel: 

fi'a + 9b, , 

we have 


C 


cd 

lab 


rjia rjia 
ac^ cd 


C ed _ jabc' rpe' 

2ab — ^ -^cd ^ 


'icd _ ^ 4^abc rj-}. 


r^cd _ • £Q 

‘-'Safe — 


cd * 


That will also apply to qg —?■ qg channel, 

q^ + 9a ^ qd + 9b, 


where we have 


(116) 

(117) 

(118) 

(119) 

( 120 ) 


pyac _ rpa rpa pyac _ Jobc rpc pyac _ • fabc rpc 

^Ibd — -^ac-^cd ^ ^2bd ~ ^ cd ^ ^3bd — cd ' 

For gg —)• gg channel, 

9a + 9b ^ 9c + 9d , 

however, it is much more complicated, 

^afe ^^pbe' ^ede' ^abc' pdc' (J°'^ -|- 

4 8 

Q 1 

p^ab _ jacc' jbdc' p^ab _ ^ £acc' fbdc' 

^5cd — J : 

- Ipcc'jbdc' ^ 

2 3 

Ci = + i5“‘‘i'“) - , 


( 121 ) 

( 122 ) 


(123) 
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where we have eight independent color bases. 

With the above color bases, we can decompose the soft factors in the matrix form. The 
associated soft gluon radiation is represented by eight gauge links. This is because all the 
initial and final state can radiate or absorb soft gluons. Therefore, we can decompose the 
soft factor, according to the following formula, 

SfAh) = f 

° r(i-,) 

x4„(6±)£««(6x)4H(0)’C»i',(0)|0) , (124) 

where IJ represent the color indices in the color-space bases constructed above. Therefore, 
Sij is the matrix element in the associated color-space bases for a particular partonic chan¬ 
nel. As mentioned above, we have four gauge links associated with two incoming partons, 
for which we follow the TMDs to adopt off-light-cone vectors v and v to construct the gauge 
links. The off-light-cone vectors are applied to regulate the light-cone singularities. For the 
two outgoing partons, we apply the off-shellness to cast the out-of-cone radiation contribu¬ 
tion to the soft factor. This regulation depends on the jet size. That is why we introduce 
the off-shellness = R\P^/Q‘^ and = R\P^/Q^ for the two hnal state jets with jet size 
Ri and i? 2 , respectively. Cjk/ and CJhi represent the corresponding color conhgurations, as 
introduced above. 

It is straightforward to calculate the soft factor at the leading Born order, Su = 

For the channels in Eqs. (1115111131) . we hnd 


^(0) 



0 

C/iCf 

2 


For the channels in Eqs. flll6llll8lll20p . we have 


(125) 


^(0) 


2CICf 0 0 

0 {C\ - A)Cf 0 


(126) 


0 0 CICf 

For the channels in Eq. (11221) . we obtain 


^(0) 


"5 00000 0 0" 
050000 0 0 
005000 0 0 
0 0 0 1 0 0 0 0 
000080 0 0 
000008 0 0 
000000 20 0 
000000 0 27 


(127) 


The hard factor H should be expanded by the same color bases. At the tree level, our 
results are consistent with those in Ref. [11]. For completeness, we list these results in the 
following. For the partonic channel of + (f ^ + (f, 


rW 


Hu 

H21 H22 


( 128 ) 
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where 


Hu = 


2Cl 


2CFf + u^ 2 Cf 

+ ^Vt ’ 

2 4 4 

22 — ^4 7 ^ 1 “ ^2 I2 /^3 • 


( 129 ) 


Similarly, the partonic channel —)■ q^ + q^ is expressed in the 2x2 matrix with 

^ _2Clt‘^ + u^ 

-nil — -- 


Hi2 = H21 = - 


2Cf 


H22 - ^ 


2 t^ + u^ 


(130) 


And for channel g^ + q^ —)■ g^ + q^, 


Hii — Hi2 — H 21 — 0 , 
2 


H22 = 


Ci 


(131) 


For channel g^ + g'^ —)■ g^ + g”, 


Hu = 


2Cl e + s2 


u - U - ‘^Cf^ + s^ , 2Cfs^ 


Hoo = 


2 + u^ 2 + r 2 

^ t2 


(132) 


For channel g^ + q^ —)■ g^ + g-^ , 


FTii — Hi2 — H 21 — 0 , 


H22 = 


2 s^ + u^ 


c'i 


(133) 


For channel g + g ^ q + q and g + g ^ g + g, the hard factor is calculated as 3 x 3 matrix, 


i3'»' = 


Hu 

Hi2 

CO 

H21 

H22 

H23 

Hs, 

H32 

CO 

CO 

28 




(134) 



where 


^11 

Hi2 

H22 

Hi3 

H23 

H33 


1 

2C\ ut 


H2^ = 


1 + r 

2C\ ut 


1 u‘^ + t‘^ 
2C\ ut 


H31 = 


H32 = 


1 


1 t^ —u^ 
2C\ ut 
1 t^ —u^ 
2C\ ut 
4 tu 


1 t — u 

1 t — u 


2Cl ut Cl s2 Cl ■ 


(135) 


For channel q + g-j-q + g, we hnd a similar 3x3 matrix with 


ffii 

H32 

H22 

H33 

H23 

H33 


1 C — 2su 

AC\Cf ’ 

1 — 2su 

^ 4C'iC'i. su ’ 

1 t^ — 2su 

AClCp ’ 

I sH + As^u + 2stu — tu^ 
2 CICf stu 

1 + 4s^m + 2stu — tu^ 

“ ~2CICf stM 

1 {2su — — Asu) 

ACICf ■ 


For channel g + g ^ g + g, however, we have 8x8 matrix, 


(136) 


'0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Hm 

H 45 

H^e 

0 

H 48 

0 

0 

0 

H 34 

H 55 

Hu 

0 

Hu 

0 

0 

0 

Hu 

Hu 

Hu 

0 

Hu 

0 

0 

0 

0 

0 

0 

0 

0 

.0 

0 

0 

^84 

Hs5 

Hse 

0 

Hs 8 


29 


( 137 ) 



where 


Hu 

Hu 

H55 

Hu 

Hm 

Hm 

Hu 

Hu 

Hu 

Hs8 


9 (s^ — tu){st — u^) 


16 

Hu = 


s'^u^ 


9 (s^ — tu){st — V?) 


32 


9 9 ’ 


9 (s^ — tu){st — u^) 


64 

Hu = 
Hu = 




9 {s — u) (s^ + SM + M^) {su — t^) 
9 {s — u) (s^ + SU + u^) (su — t^) 


64 s'^tu^ 

9 (s^t — 2su^ + t^u) (2s‘^u — st^ — tu^) 


Hu = 


Hs5 = 


Hse = 


QAsH'^u'^ 

3 (s^ — tu){st — u^) 


s'^u^ 


16 

3 (s^ — tu)(st — u^) 

^ ’ 

3 (s — u) (s^ + SU + u^) (su — t^) 


32 sHu^ 

1 (s^ — tu)(st — u^) 


16 


s'^u"^ 


(138) 


The above hard factors are normalized to reproduce the leading order differential cross 
sections in Sec. II, 


,{0) _ c.l2^34(0) ^12^34(0) OQ^ 

^12^-34 ~ ^IJ Hjj , (139j 

where sum over JJ is understood. Because the leading order soft factor is diagonalized, it 
is simple to verify the above equation. 


C. Soft Factor and the Anomalous Dimension at One-loop order 

At one-loop order, we will expand the soft factor dehnition with perturbative corrections. 
The gluon radiations between all gauge links will contribute. For convenience we separate 
out the common kinematic integrals from the color factors for each soft factor calculation, 

S“> = W„I.i . (140) 

ij 


where ij label the associated gauge links: 12 for the incoming two partons and 34 for the 
outgoing partons. In the above equation, lij represent the kinematic integrals for the soft 
gluon radiation between i and j gauge links, whereas the factor Wij represent the associated 
color factor in the matrix form for particular partonic channel. 

We show in Fig. ([7]) a few representative real gluon emission diagrams, labelled by the 
relevant kinematic integral Jjj, which contribute to the soft function at the next-to-leading 
order (NLO) in as- 
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FIG. 7. Representative diagrams contributing to the soft function at the NLO. 


For the -^34 diagram, with gluon connecting 3 and 4 gauge links, we find the following 
results: 


(sin0) 

7^r(i-e) 

r(i-e) 


as J_ 
27i‘^ 


in the transverse momentum space, where k represents the radiated gluon momentum. We 
have kept the e term in the above equation, because they will contribute to a finite term 
when Fourier transforming to fe^-space. In the calculations, tij and fij carry the momentum 
direction of two outgoing jets, and therefore, the scalar product nj-fij enters in the final result 
which is proportional to Q‘^/P‘1. That is the cause of having the logarithmic dependence 
on this factor in our final result. Again, the offshellness of the two outgoing gauge links 
are introduced to cast the contribution that is only associated with out-of-cone radiation. 
Gluon radiation within the jet cone will not induce a nonzero transverse momentum of the 
dijet and its contribution needs to be subtracted out. 

Similarly, for gluon connection between all other gauge links needed in the final calcula- 


In 




RlRjP^ ) 


+ e 


1 , 2 ! 1 , 2 ! ^ 5, s 

- In^ — 7 T H— In^ —-I-4 In — In- 

2 Rf 2 Rl 3 -t -u 


dk^dk 


2nj ■ Uj 


{k ■ nj){k ■ rij 




(141) 
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tions, we find 
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(143) 


(144) 


(145) 


(146) 


(147) 


Here, we did not list either /as or 144 , because they will not contribute to the soft factor. 
For the contribution from virtual diagrams, all the integrals are proportional to 

2 -^ . ( 148 ) 

^IR ^UV 

After introducing the UV counter-terms which are proportional to 

-+ 7 ^ + In ^ , (149) 

€uv 47r 

only IR divergence will be left and cancel the singularity in real diagrams, after we have 
performed Fourier transformation from the transverse momentum space to space. 
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The final expressions in the 6_L-space can be written as 
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(151) 

(152) 

(153) 

(154) 

(155) 
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The color factors can be calculated diagram by diagram as well, and for the the channel of 
g + g —)■ g + g, they are: 


VTli = 1^22 = Cf 


W 34 = W,2 = 

fTi3 = 14^24 = 

1^23 = Wu = 


n CaCf 
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For the channel of g + g —)■ g + g, 
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For the channel g + g ^ q + we have 
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For the channel q + q ^ g + g, 
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For the channel q + g ^ q + g, the color matrixes are, 
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For the channel g + g ^ g + g, they are 
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we 

will be able to calculate 


(181) 


(182) 


(183) 


soft factor, 




(184) 


for all partonic channels, where nm run from 1 to 4. It is important to note that after 
summing the contribution in every diagram, the jet cone size R 12 and p dependence will be 
canceled by each other in all the non-diagonal elements in and their color factors are 
probational to For example, for the contribution from q + q ^ q + q channel, we hnd 
the soft factor can be written as 
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in the transverse momentum space, where, for convenience, we have dehned 
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In the above equation, we have also introduced a short notation for an additional 

term, which will be dehned later, together with those for all other channels. 

For g + g —)■ g + g, we have 
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For g + g —)■ g + g, we have 
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For g + g ^ g + g, we have 
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And finally, for g + g —)■ g + g, we obtain 
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In the above equations, the S matrixes for all different channels are dehned as 
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Similarly, we obtain the soft factors in the 6 _L-space, 
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where Inmip) have been calculated in Eqs. fll501ll56p . For channel g + g g + g, we find 
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For g + g —)■ g + g, we find 
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For g + g ^ g + g, we find 
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For g + g ^ g + g, we find 
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For g + g —)■ g + g, we find 
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For g + g —)■ g + g, we find 
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We would like to emphasize again that the nontrivial results of the soft factor calculations 
in the above. In particular, the p and 2 -dependence only appear in terms proportional 
to This is consistent with the factorization we have argued in the beginning of this 

section. We will show explicitly how the factorization works in the following sub-section. 


D. Factorization at One-loop Order 

In this subsection, we will apply the one-loop results for the soft factor and TMDs calcu¬ 
lated in the last subsection to verify the factorization formalism we have proposed. We will 
show the factorization in both transverse momentum space and Fourier conjugate 6_L-space. 
At his order, the TMDs can be calculated, and have been listed in Sec. VIB. The soft factors 
have also been calculated above in Sec. VIE. 


1. Transverse Momentum Space Factorization 


We hnd from Eq. (11021) that the dijet differential cross section at a nonzero transverse 
momentum q± receives contributions from the two incoming parton distributions and the 
soft factor. The hard factor, for a given 2 —)■ 2 process, does not depend on as expected. 
Therefore, we can write down the expansion of the hnite q± contribution from the TMD 
factorization formula at one-loop order: 
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where we have expanded the TMDs and the soft factor at one-loop order. We have calculated 
all these factors in previous subsections. Substituting these results into the above equation, 
we can derive the hnite q± expression for dijet production via each channel, and reproduce 
the results shown in Sec. IV. Especially, the p-dependence from both the TMDs and the soft 
factors is cancelled out by each other, which is a nontrivial illustration of the factorization 
formula. 


2. Factorization in b±-space 


Similarly, the factorization can be demonstrated in the 6_L-space, for which we carry out 
the calculations for W{b±) in Eq. fllOSp . 
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With the expressions calculated in previous sub-sections, we will be able to reproduce the 
logarithmic terms in Sec. IV. In principle, we should be able to extract the hard factors at 
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one-loop order as well. However, to do that, we need the expressions of the virtual diagrams 
contributions in the color-space constructed in Sec. VIC. We hope these calculations can 
be performed in the future, from which we shall obtain the hard factors at one-loop order. 
However, to show the factorization, we do not need the hard factors at the one-loop order. 


E. Resummation 


According to the dehnition of the W function, cf. Eq. fllU3p . when we choose factorization 
scale /i = there will exist two classes of scales in TMDs and soft factor, one is which 
comes from the Fourier transformation of transverse momentum hence it is a small scale, 
the other class contains scales of and where CiCl = ^-re large scales. Our 

one-loop results for the TMDs and soft factor show that there are large double and single 
logarithms of InQ^fe^ in them. (Here, we use h± and b exchangeably.) These large logarithms 
appear at every order of perturbative expansion as the form of a” ln^"'Q^6^ or a” ln”Q^6^, 
when In Q‘^b‘^ > l/a*, the convergence of the conventional perturbative expansion is impaired 
and no longer gives a correct prediction. In order to make a reliable calculation for the W 
function, these large logs have to be resummed. Based on the factorization theorem, we 
could hnd all the relevant scales present in each factorized factor, which can be evolved 
individually and satishes a renormalisation group equation(RGE). For example, for the 
TMD f{x,bX,P,p), the relevant Collins-Soper evolution equation reads as: 

b, C, /i, p) = {K{b, p) + G(C, p)) f{x, b, C, /i, p) , (206) 

According to the one-loop result, we can get 

K{b, fi) + G(C, fi) = -^ In . (207) 

71 4 

where Cj = Ca for gluon distribution and Cj = Cp for quark distribution. K{b,p) satishes: 
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The renormalization group equation for K{b,p) is governed by the associated anomalous 
dimension jk, 

dK{b,p) dG{C,p) 2asCi 

d\np d\np tt 

Then we can solve the above renormalization group equation. 



K{b, p) + G(C, p) = K{b, Qo) - r %1 k{p') - ^ In , (210) 

jQo P 71 p 

where in the end we will choose Qo = b^/b to resum the large logarithms. In addition, we 
can evolve the factorization scale in the TMD distribution as well. 
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with the anomalous dimension 
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Here, for gluon distribution we have Bj = (3 q, and for quark distribution we have Bj = 

By solving the above evolution equations, we resum the large logarithms associated with 
the TMDs, 
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where we have chosen = pQ^. All the large logarithms have been resummed into the 
Sudakov form factors. In addition, applying the one-loop results in Eqs. fllOSIIlOSp . we find 
that the large logarithms (In associated with the integrated parton distributions can 
be resummed by choosing the relevant scale as p = b^/b. By doing that, the hrst factor in 
the above equation can be replaced with 


*>. A. p p) ^ /(a:, 1* = (1 + a, ■ ■ ■) . 


(214) 


where the right hand side is the integrated parton distributions at the scale oi p = bo/b, 
and we have also neglect all the constant terms, such as p-dependent terms. These terms 
are beyond the accuracy of NLL we are considering in this paper. 

For the soft factor, it satishes the evolution equation of 
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Here T'® is an identity matrix, for di-gluon initial states A = Ca, for di-quark initial states 
A = Cp and for quark-gluon initial states A = {Cp + Ca)/2,. For quark jet Dj = Cp, while 
for gluon jet Dj = Ca- According to the one-loop correction of the 7 ^ factor for each 
channel can be directly read out. For example, for the production channel q + q ^ q + q, 
we hnd 
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For the channel q + q ^ q + q, 
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For the channels g + g ^ q + q and q + q ^ g + g-, 
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By solving the evolntion eqnation, we obtain 
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, ( 222 ) 


where the hrst factor conies from the contribntion of the F^ term in Eq. fl216p . We note 
that since there is only one scale present in S{b,qi = bo/b), it does not contain any large 
logarithm. 

Snbstitnting the above solntions, Eqs. fl213p and fl222p . into the factorization formula, 
Eq. fllOSp . we obtain the hnal resummation results for W. In particular, since we have 
set the factorization scale g, = Q, there will be no large logarithms in the hard factor H. 
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Therefore, all the large logarithms have been resummed into the Sudakov form factors, and 
we have 


W (xi, X2, b) = Xi fa{xi, bo/b±)x2 fb{x2, bo/b±)e •^sud(Qh 6 ±) 
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Jbo/b^ Jbo/b^ 


X Tr 


with 


^sud(Q^&±) = 






hl/bl 


In I 71 +B + 


which are exact the results we showed in the Introduction section. 


(223) 


(224) 


F. Contributions from the Non-global Logarithms 

In Refs. [3,1^, the so-called non-global logarithms were discussed for the dijet correlation 
in hadronic collisions. They further take an example of dijet production in DIS processes, 
and estimate their contributions [^. This non-global logarithm comes from the kinematics 
of two gluon radiations, where one is within the jet and another soft gluon outside the jet. 
Since it happens at C>(a^) order, the one-loop calculations in this paper do not encounter 
this non-global logarithm □. Numerically, the non-global logarithms are negligible. This is 
because it starts at and in the kinematics of low imbalance transverse momentum, 

the resummation is overwhelmingly dominated by leading double logarithms contributions. 
In the following calculations, we will not consider their contributions when we compare to 
the experimental data. 


VII. PHENOMENOLOGY OF DUET CORRELATIONS AT TEVATRON AND 
LHC 

In this section, we will apply our resummation formula to the dijet production from 
collider experiments, including Tevatron and the LHC. In these experiments, the leading jet 
energy is large, and we expect that the resummation is dominated by the perturbative form 
factors. That means that our predictions are not sensitive to the non-perturbative form 
factors. 


A. Non-perturbative Form Factors in the Resummation 

To apply the resummation formula for phenomenological applications, we follow the 6*- 
prescription to introduce the non-perturbative form factors j^, he., 

-)■ 6* = b/ a/I + ! (225) 

^ It seems possible that this contribution might belong to the soft factor in our factorization formula of 
Eqs. (I102I103P at two-loop order. At this order, the soft factor contribution will have similar kinematics 
as described in Q for the non-global logarithm, where one gluon is within the jet (the Wilson line in the 
soft factor definition) and one gluon is soft and outside the jet. Needless to say that this has to be verified 
by an explicit calculations with two gluon radiations. We plan to come back to this issue in the future. 
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in the 6-space cross section contribution W{h_i) with b^ax a parameter which will be set as 
bmax = O.SGeV”^. By doing that, it is guaranteed that 6* is always in the perturbative 
region. Therefore, W{b±) is replaced by 

W{xi,X2,b) —)■ W{xi, . (226) 

The non-perturbative form factors follow the parameterizations in Ref. jl^. Since we have 
quark and gluon from the initial state, we decompose the non-perturbative form factor into 
the quark and gluon contributions, 

b) = 6) + i>) . (227) 


for a partonic channel ab —)■ cd. In the right hand side of the above equation, Snp depends 
on the flavor of the incoming partons, 

S'np = ^ 1" (^) + 9s’’ lii(10ii)6’ , 

^NP = 4"''' + 4" (^) • (22®) 

with the following parameters: 

g[^^ = 0.21, (7^"^ = 0.68, g^^^ = -0.29, 

g[^'> = 0m, # = 0.87, # = -0.17. (229) 


These parameters are htted to the hard scattering processes in the relevant qq and gg 
processes in Ref. [1^. In our calculations, we assume that these non-perturbative form 
factors apply to the dijet production processes as well. This is an approximation. However, 
we would like to emphasize that because the jet energy is so large that our hnal results are 
not sensitive to the non-perturbative form factors at We have also checked several recent 
proposals for the non-perturbative form factors [ 


261], and found that all of them predict 


almost the same distribution for dijet production in the following phenomenological studies. 


B. Dijet Correlations at the Tevatron 


With all the ingredients calculated above, we now compare our resummation results to 
the experimental data from the Tevatron. In experiments, the normalized differential cross 
sections are measured, 

1 dcTdijet 



where 0dijet is the azimuthal angle between the leading two jets. The leading jets are in 
the separate transverse momentum bins, with the second leading jet transverse momentum 
Pt > 40 GeV. The events are selected in the mid-rapidity, \yjet\ < 0.5. 

In Fig. [HI we plot the comparisons between our resummation result and the experimental 
data from DO collaboration at Tevatron. For completeness, we also show the prediction of 
a hxed order calculation at the NLO , which includes both one-loop 2 —)■ 3 and tree level 
2 —)■ 4 contributions. To compare with the normalized differential cross section, we have 
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FIG. 8. The comparisons between the theory calculations and the experimental data from the DO 
collaboration at the Tevatron. 

normalized the result of our resummation calculation in 0dijet distribution by the LO dijet 
cross section, with both jets in the specified pt and y bins of the data point. Namely, the 
normalization factor Udijet is taken to be the NLO dijet cross section in the NLO predic¬ 
tion, and the LO dijet cross section in the resummation prediction. This is because in our 
resummation calculation, cf. Eq. fl223p . the hard factor H is only kept at the LO in this 
calculation. From these plots, we can clearly see that the resummation results agree well 
with the experimental data around the back-to-back correlation region of 0 ~ tt. For smaller 
value of (f) (away from the back-to-back configuration), the resummation calculations match 
to the fixed order results at NLO [^, which has also been separately shown in F^. [HI We 
note that a full NLO calculation cannot describe experimental data for 0 ~ tt [ij, where 
the fixed order calculation becomes divergent. Our resummation calculation, after being 
matched with the NLO result clearly improves the theory prediction and can describe the 
experimental data in a wider kinematic region. This demonstrates the importance of all 
order resummation in perturbative calculations for these type of hard QCD processes. 


C. Dijet Correlations at the LHC 

Dijet production processes are among the hrst few measurements of pp collisions at the 
LHC. Both CMS and ATLAS have reported the experimental results on the azimuthal 
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CMS VS = 7TeV 

300 < P™’' < 1000 (GeV) (XI000) 
200 < Pj®’' < 300 (GeV) (XI000) 
140 < Pj®’' < 200 (GeV) (XI00) 
110< Pj®’‘< 140 (GeV) (XI0) 

80 <Pt <110(GeV) 

1= anti-k, R=0.5 > 30GeV |y|<1.1 



FIG. 9. The comparisons between the resummation results and the experimental data from the 
CMS collaboration at the LHC. The kinematics of each bins are specified according to the experi¬ 
ment [3]. 


angular correlations of dijet productions, as done by the DO Collaboration at the Tevatron. 

In Fig. [9l we compare our resummation results with the experimental data from the 
CMS collaboration at the LHC. Similar to the DO measurements, the dijet measurements 
are presented in several kinematic bins, with the leading jet transverse momentum labelled 
by in the figure. The second jet transverse momentum is chosen to be larger than 

30 GeV. Both jets are in the mid-rapidity region, \yjet\ < 1-1. Anti-fc* jet algorithm with 
jet size R = 0.5 was used in the data analysis. We have also applied this algorithm in 
our calculations. In this figure, we limit the comparisons in the back-to-back correlation 
region, where we find perfect agreements between the resummation calculations and the 
experimental data over all transverse momentum bins. Similar to that in Fig. [HI away from 
the back-to-back region, the resummation calculations will match to the fixed order results. 

In Fig. Iini we compare to the measurements form the ATLAS collaboration j^. In this 
experiment, the same anii-kt algorithm has been used, however, with jet size R = 0.6. The 
two jets are selected from the mid-rapidity region {\yjet\ < 0.8) with minimum transverse 
momentum of 100 GeV. The data sets are chosen according to different leading jet transverse 
momentum region as indicated in the figure. From this figure, we can see that the agreements 
between the resummation results and the experimental data are very well around the back- 
to-back correlation regions, except in the lowest bin. The apparent poor agreement 

between the resummation prediction and the ATLAS data in the lowest bin (between 
110 GeV and 160 GeV) is caused by the stronger kinematic cut made on the second jet Pt, 
which is required to be above 100 GeV at ATLAS and 30 GeV at GMS. With a much tighter 
cut on this second jet Pt, the phase space for multiple soft gluon emission is limited so that 
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FIG. 10. The comparisons between the theory calculations and the experimental data from the 
ATLAS Collaboration at the 7 TeV LHC. 

our resummation calculation (which allows all possible soft gluon radiation) becomes less 
reliable in this case. We note that in the lowest bin, the cross section is dominated 

by around 110 GeV which is close to the 100 GeV cut on the second jet made by the 
ATLAS. 

VIII. SUMMARY AND DISCUSSIONS 

In summary, in this paper, we have investigated all order soft gluon resummation in dijet 
production processes in hadronic collisions. The procedure and methodology follow the orig¬ 
inal GSS resummation for massive neutral particle production. Because the hnal state jets 
carry color, the resummation formulas have to be modihed to include the effect of soft gluon 
radiation associated with the hnal state jets. In the derivation, we calculated the complete 
one-loop contributions from soft and collinear gluon radiations in all partonic channels in 
dijet productions. The soft divergences are shown to be cancelled out completely between 
real and virtual graphs, which provides an important check to our one-loop calculations. In 
order to derive an analytic expression with the jet cone size dependence to demonstrate the 
cancellations in the hnal results and to derive the resummation formula, we apply the narrow 
jet approximation in our calculations. We have also implement the (anti-/Ci) jet algorithm 
to separate out the out-of-jet cone radiation from the gluon radiation inside the cone jet. 
Hence, the hnal results of our calculation depend on the jet algorithm and the jet size, which 
can then be directly compared to experimental data. As an important cross check, we have 
compared our derivation of soft and collinear gluon radiation contribution to the hxed order 
calculations at this order, and the numeric comparisons show that they agree very well in 
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the kinematics of back-to-back correlation regions. In this region, because the predictions of 
hxed order calculations are divergent, we have to take into account all order resummation 
effects. 

We have compared our resummation results to the experimental data from the measure¬ 
ments at the Tevatron and LHC. All these comparisons demonstrate that the resummation 
results are crucial to improve the theory descriptions of the experimental data around the 
azimuthal back-to-back kinematic region. The combination of the NLO perturbative cal¬ 
culations (including both one-loop 2 —)■ 3 and tree level 2 —>■ 4 contributions) and our 
resummation results provide the most adequate theory descriptions to these experimental 
data. 

Our calculations are the hrst systematic derivations of the TMD resummation for dijet 
production in hadronic collisions at the NLL order. The results have been cross checked 
through various perspectives, and they are consistent within the theoretical framework our 
calculations are built on. These cross checks are nontrivial supports for the factorization 
arguments used in our derivations. A number of extensions can be performed along this 
direction. For example, we shall be able to calculate the soft gluon resummation effects 
in the vector boson (or Higgs boson) plus a high Pt jet production at the colliders where 
the total transverse momentum of the boson and the jet is much smaller than the invariant 
mass of the hnal state particles jl^. These processes are important channels to study the 
Standard Model physics at the LHC. 

Finally, we would like to comment on the applications of our results to the dijet production 
with large rapidity separation. This particular kinematics is very interesting to study the 

S D resummation physics. It has long been realized that the so-called BFKL resummation 
will be important in this kinematics, which is referred as the Mueller-Navelet dijet 
production [29|. In particular, recently, the CMS collaboration at the LHC has measured 
the dijet azimuthal correlation with large rapidity separation between the jets, which has 
been interpreted as the BFKL resummation effects |30|. However, in this calculation, only 
BFKL-type resummation has been taken into account. We would like to argue that there 
should be Sudakov resummation as well. Theoretically, how to resum large logarithms from 
both types of physics effects is an important question. We will not address it in this paper. 
Instead, we will discuss below the physics of our resummation formula when considered in 
this kinematics. 

When the two jets are produced with large rapidity separation, we are in a special kine¬ 
matic region, where the physics is dominated by f-channel diagrams. Therefore, we can 
apply the following kinematic approximations, s ~ —u S> —t, which also implies that 
= tu/s ~ —t. More importantly, all the partonic channels with t-channel gluon ex¬ 
change will be the most important contributions. This is because they all have terms which 
are proportional to This includes the following channels: qq' —)■ qq', qg —)■ qg, and 

99 gg. In addition, by applying the above approximation, we find out that the anomalous 
dimensions for the associated soft factors derived in the last section become diagonalized. 
The direct consequence is that we can simplify the final resummation formula, by absorb¬ 
ing the soft factor anomalous dimension of Eq. fl223p into the overall Sudakov perturbative 
form factor of Eq. fl224p . This much simplihed result, as compared to that presented in the 
last section, may indicate a consistent resummation formula for the Mueller-Navelet dijet 
production. The remaining task is to develop a consistent theoretical framework to include 
both physics effects induced by the BFKL and Sudakov resummation dynamics [3l|. 
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